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ABSTRACT 


Heat  and  mass  transfer  studies  for  flow  over  a  cylinder  having  a 
longitudinal-strip  source  of  heat  or  mass,  are  presented.  A  theory 
for  the  two  asymptotic  cases  of  P”  oo  and  Pr  ~*0  is  developed,  giv'ing 
the  temperature  profile  and  the  transfer  rate  for  such  a  strip.  The 
theory  includes  the  effect  of  the  conductivity  of  the  heating  strip  mater¬ 
ial.  It  is  shown  that  if  the  conductivity  of  the  strip  material  is  v’ery 
high,  then  the  transfer  rate  is  a  measure  of  local  wall  shear-stress 
provided  the  Prandtl  number  of  the  convecting  fluid  is  high. 

Experimental  studies  described  include  electrochemical  mass  - 
transfer  studies  carried  out  on  a  rotating  cylindrical  electrode  carrier; 
and  heat  transfer  studies  carried  out  in  an  oil  tunnel.  Fair  correlation 
between  theory  and  experiment  w'as  found. 
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1 .  Introduction 

The  purpose  of  the  present  report  is  the  presentation  of  both  theoretical  and 
experimental  results  for  ths  transfer  of  heat  and  mass  from  or  to  a  long  circular 
cylinder  immersed  in  a  two-dimensional  flow  of  an  incompressible,  isotropic, 
Newtonian  fluid.  The  flow  is  of  laminar  boundary  layer  type,  with  a  step  type 
heat  (or  mass)  flux  source  on  the  cylinder  envelope. 

The  subject  of  heat  transfer  in  external  flow  has  been  studied  previously  by  a 
large  number  of  investigators,  both  th'^oretically  experimentally  ^8,  *. 

In  particular,  asymptotic  solutions  for  the  cases  of  e.xtremely  large  and  extremel'/ 
small  Prandii  or  Schmidt  moduli,  and  for  many  kinds  of  temperature  boundary- 
conditions  at  the  cylinder  surface  ha\  e  been  known  for  some  time  [,’^3  .  Now,  hese 
analytical  solutions  depend  upon  the  detailed  knowledge  of  the  e"tcrrial  flow  around 
the  cylinder  as  a  function  of  coordinate.  This  information  is  ob‘.-..n ..ble  cither 
from  the  analysis  of  an  assumed  idealized  flow  field  (see  c.f.  ^loj),  at  'east  up 
to  that  value  of  the  coordinate  in  the  flow  direction  where  boundary  layer  separation 
takes  place;  or  else  it  has  to  be  obtained  from  a  semi-empirical  approach  which 
takes  into  account  actual  test  conditions  (see  c.f.  ',13,  14]  ).  However,  beyond  the 
point  of  boundary  lay^er  separation  little  is  known  about  wall  shear-stress  values 
and  local  heat  and  mass  transfer  rates. 

There  has  recently'-  been  some  interest  in  the  reir.vestigation  of  two  dimensional, 
incompressible,  laminar,  steady  and  separated  flow  behind  bluff  objects  Cl4.  ’3]  . 
This  problem  is  both  of  obvious  practical  significance  in  the  estimation  of  the  drag 
of  such  bodie-',  and  of  fundamental  theoretical  importance  in  the  investigation  of 
the  stability  of  two  dimensional  flows  and  the  viniqucncss  of  the  flow  field  at 
(comparatively)  high  values  of  the  Reynolds  number.  The  choice  of  suitable  flow 
models  tc  represent  such  separated  flows  must  perforce  depend  on  at  least  some 
rudimentary  knowledge  of  the  flow  conditions  inside  the  region  of  separation  (the 
wake  bubble). 


*For  a  summary  of  previous  work  see  c.  f.  LlOi  1  ll  and  [l2  3* 
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The  present  report  will  describe  a  mass  and  beat  transfer  experiment  in  which 
the  flux  •nurce  (or  sink)  is  a  riarrow  strip  bounded  by  tv/c  generating  lines  on  the 
cylinder  envelope.  It  will  be  shown  that,  provided  this  strip  is  ?u.ii';iently  uarrow, 
the  transfer  rate  from  its  surface  depends  only  on  the  local  wall  shear  stress. 

Thus  the  measurement  of  the  transfer  rate  may  be  translated  immediately  into  a 
local  wall  shear  stress,  information  probably  not  obtainable  from  any  other  simple 
experiment. 

The  design  of  the  experitnental  apparatus  will  be  desci'ibed,  and  expe.Jiime’itai 
results  will  be  compared  with  theoretical  predictions. 


2.  Theory 


It  will  be  found  useful  to  review  here  briefly  Lighthill's  [  3  ]  asymptotic  solution 
for  the  heat  (or  mass)  transfer  to  laminar  boundary  layer  flows  at  very  high  and 
very  low  values  of  the  Prandtl  modulus  (respectively  the  Schmidt  modulus  in  the 
case  of  mass  transfer). 

We  shall  introduce  the  boundary  layer  approximations  into  the  dimensionless 
equations  of  momentum,  continu  ity  and  energy  for  the  flow  of  an  incompressible, 
non-dissipative,  isotropic,  Newtonian,  constant  physical  property  fluid  to  obtain. 


8u  ^  ^  9*u 

9x  9Y  ~  ‘  '9x  9Y* 


^  .  9V 
9x  '  9V 


0 


^  V  ~  —  — 

9x  ^  ^  9Y  ■  Pr  8Y*" 


with  boundary  conditions, 


X  <  x^;  e  -  0; 

x  -  0,  Y  =  0;  u  ^  V  =  0; 

Y  =  03;  9  =  0;  u  =  U(x); 

X  -  x*;Y  =  0;  G  =  9(x):'  ' 

In  what  follows  interfacial  velocities  w'ill  be  neglected  so  that  the  equations 
above  will  be  applicable  to  both  the  transfer  of  heat  and  mass.  For  the  latter  case 
the  Sherwood  number  is  to  be  read  whenev'er  the  Nusselt  number  appears,  and  the 
Schmidt  number  for  the  Prandtl  number. 

2.  1  Large  Prandtl  or  Schmid  moduli 

For  Pr  the  thermal  boundary  layer  may  be  assumed  to  be  very  thin  in 
comparison  with  the  momentum  boundary  layer.  u(x,  Y),  being  an  analytic  function 
within  the  flow  field,  may  thus  bt  expanded  in  a  Taylor  series  arovind  its  value  on 


(*)  see  list  for  explanation  of  symbols 
(:,■;*)  respectively  90/9y  prescribed. 
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the  boundary.  Following  Leveque  [,163  only  the  first  term  of  the  expansion  will 
be  retained  (see  Appendix  I).  Then  a  solution  of  equations  (Z.  1)  satisfying  (Z.  Z) 
is  found  to  be  (see  c.f.  [7'J  ), 


9 


^  /  -'i'  H 

■  ^ 

391  _  XT  _  Pr*^ 

dy  y=o  ■  ■  9  "  r  C  j) 


d  9s (x) 


\ 

! 


i.  ...(2.3) 


✓ 


Here  x*  is  the  dimensionless  abscissa  of  the  first  point  of  non-zero  9^  ■ 

Several  sv^bcases  representing  idealized  situations  may  now  be  considered: 

2.  1-1  Interface  temperature  (or  concentration  in  the  case  of  mass  transfer)  zero 
for  x<x*,  and  equal  to  a  non- zero  constant  for  x  -  .x=:^ 

For  this  case  the  second  of  equ'dions  (2.  3)  reduces  to, 


Nu  = 


Re^  Pr 

9‘<j  r{%) 


k:  vAj)  =15  ] ' 

To  get  the  average  value  of  the  Nusselt  modulus  in  two  dimensional  flow  over  a 
heating  surface  extending  from  x-  to  x,  we  may  integrate  equation  (2.  4): 

3''^  -  V  n  d^] 


(2.4) 


Nu  = 


Re 


.  .  .  (2.  5) 


If  the  heated  surface  is  relatively  narrow  in  the  x- direction,  then  may  be 

assumed  to  be  essentially  constant  from  x*  to  x;  the  equation  above  then  simplifies  to  - 


where  ^  is  some  average  value  of  : 


)  = 


X  -  X^ 


...(2.6) 


...(2.6-a) 


It  is  apparent  from  equatioil  (Z.  6)  that  the  measurement  of  the  surface  temperature  o 
a  strip  of  suitable  dimensions  and  the  calculation  of  the  average  dimensionlt  ss  flux 
will  furnish  directly  an  amplified  value  of  the  local  shear  stress  at  the  wall,  quite 
independently  of  previous  boundary  layer  history. 
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2.1-2  Zero  interface  heat  flux  for  x<x^,  non-zero  interface  flux  for 

We  shall  investigate  the  case  in  which  the  surface  flux  changes  in  steps  of  finite 
size:  cases  of  continuous  flux  variation  are  easily  derived  from  these  results  by 
superposition  and  passage  to  the  limit.  Define, 


Ai(x,)  =  q."  - 


,  L  Pr  Re' 


—  (i  =  0.  1,2 . n) 


.  .  . (2.  7) 


1  \  4  ! 


where  q/'  is  a  surface  heat  density  increment;  1  is  a  characteristic  length  and 
x^  is  the  coordinate  location  considered.  Then, 

0g(Xj)  =yVo(x.,)  +  (^7V,(x«  )  j  5P{x,,  -  x"0  + 

=  ...(2.8) 


X 


Here  the 


y\_  .  ar 


e  the  step  magnitudes,  occurring  at  abscissae  x^'^  and 


^  I  O' A- 


.  . .(2.  9) 


.  . . (2.  JO) 


Specializing  for  a  single  step,  we  get  Shah's  expression  ([l7],  p.  60): 
9  -  q  M  A  / _ dx _ 

Ss  (x,)  -  q;  K  30 r(-:^  J. 


Also, 


Nu  = 


0gK 


-  3  '  r(h)  Re  f  ^  — _ — _ 

1  u;pr  Re 


(2.  11) 


(2.  12) 


for  the  single  step-type  flux  variation  at  the  wall. 

The  equivalent  mean  surface  temperature  and  the  mean  Nussclt  number  ov  -r  a 
strip  of  small  finite  width  in  the  direction  of  flow,  bounded  by'  two  parallel  cylinder 
generating  lines,  may  now  be  calculated.  Proceeding  in  a  manner  analogous  to 
§  2.1-1  and  assuming  constant  over  the  strip  width  one  obtains, 


6 


and 


Nu 


...(2.  13) 

...(2.14) 


2.1-3  Cases  departing  from  the  idealized  situations  of  piecewise  constant  interface 
temperature  (respectively  concentration)  or  flux 


In  actual  practice  none  of  the  idealized  situations  analyzed  in  the  preceding 
paragraphs  is  ever  attained.  Usually  an  interface  wall  containing  distributed  heat 
sources,  cr  their  equivalent,  over  part  of  its  length  is  convection  cooled  by  some 
fluid.  In  general  the  distribution  of  these  heat  sources  v/ill  not  correspond  to  the 
rates  of  convection  to  the  cooling  fluid:  indeed,  only  one  unique  distribution  could 
be  thus  compatible.  Consequently  some  circumferential  conduction  (i.e.  conduction 
parallel  to  the  direction  <  ‘  fluid  flow)  will  ensue.  If  the  resistance  of  the  wall 
material  containing  the  heat  sources  to  the  conduction  of  heat  is  slight,  then  the 
inte-face  may  be  expected  to  approach  an  isothermal  surface.  If  however  that 
resistance  should  be  v'ery  large,  then  the  interface  flux  rate  will  correspond  to  the 
local  wall  flux  density,  and  we  approximate  case  2.  1-2.  By  small  and  large  resis¬ 
tances  to  the  conduction  of  heat  within  the  wall  m  iterial  are  here  meant  small  and 
large  values  of  the  Biot  modulus  respectively  (see  Appendix  II). 

It  will  be  shown  ir  Appendix  11  that  the  t.emperature  distribution  for  the  case  of 
an  almost  isothermal  wall  is  easily  derived  through  a  regular  perturbation  scheme, 
while  that  of  the  almost  constant  flu.x  case  may  be  solved  through  a  singular  per¬ 
turbation  method^  at  least  in  principle. 


2.2  Small  Prandtl  or  Schm.idt  moduli 


For  Pr  «  1  the  momentum  boundary  layer  will  be  rather  thin  compared  to  the 
thermal  boundary  layer.  Thus  the  momentum  boundary  layer  may  be  disregarded 
entirely  in  the  solution  of  the  energy  equation  and  the  velocity  will  be  a  function  of 
the  lengthwise  coordinate  only.  Appropriate  asymptotic  solutions  may  .hen  easily 
be  obtained. 
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2.2-1  Interface  temperature  zero  for  x<x*  and  equal  to  a  non-zero  constant  for 
xtx* 


Nu  = 


(U(x)  )' 


r^U(x)dx  1 


1 


(2.  15) 


Whence  the  average  value  of  Nu  over  a  strip  width  x^  -  x*  is  obtained, 


^  ‘‘’‘J  ' 


(2.  16) 


2.2-2  Zero  interface  flux  for  x  ^  x*;  non-zero,  constant  flux  for  x  ^  x* 

The  methods  described  in  the  previous  paragraphs  may  again  be  used.  Thus, 


NU  =  r(l) 


dx 


CvCx,)-V(>^)]^ 


.  . .(2.  17) 


with  - 


-X 

V(x)  =  J  U(^)  d^ 


We  may  again  assume  a  mean,  constant  outer-flow  velocity  to  be  operative  over 
the  strip  width;  then, 


Nu  = 


UPS' 


and, 


9 


_Jl 

K  ■  r(i) 


...(2. 18) 


...(2. IQ) 


The  further  extension  to  a  prescribed,  stepwire  variation  in  wall  flux  rates  is 
carried  out  without  difficulty,  and  will  therefore  not  be  elaborated  upon  here. 

3.  Experimental:  Mass  Transfer  Studies 

3.  1  Introduction 

An  initial  survey  of  the  relevant  literature  showed  that  the  little  work  which  had 
been  done  on  local  mass  transfer  rates  from  cylinders  in  forced  convection  at  low 
Reynolds  numbers  (Schnautz  1958,  Grassmann  et  al.  1961)  was  concerned  with 
constant  interface  concentration  over  the  entire  cylinder  envelope.  It  was  also 
found  tl.at  of  the  various  techniques  available  for  mass  transfer  studies,  an 
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electrochemical  methods  appeared  to  be  t’le  most  suitable.  According  to  previous 
studies  (  Cl9lto  j^2l3  )  such  systems  should  have  been  capable  of  yielding  results  of 
great  accuracy  and  good  reproducibility  with  a  comparatively  simple  experimental 
arrangement  and  without  elaborate  instrumentation.  Consequently  this  was  the  first 
experimental  system  chosen  for  investigation. 

3.  2  Electrochemical  details 

By  choosing  a  suitable  electrochemical  system,  in  which  electron  transfer  at  a 
test  electrode  surface  is  proceeding  at  a  very  high  rate,  it  is  possible  to  obtain  a 
range  of  E.  M.  F.  over  which  the  transfer  rate,  i.e.  the  current  measured^is  inde¬ 
pendent  of  the  driving  force  (the  E.  M.  F.  )  applied.  The  transfer  rate  will  depend 
only  on  the  characteristics  of  the  chemicals  involved,  of  the,  surface  at  which 
transfer  takes  place  and  on  the  rate  of  supply  of  ions  through  the  mass  transfer 
boundary  layer,  that  is  on  the  convection  rate.  The  current  measured  in  this  range 
is  known  as  the  limiting  current. 

Electrochemical  systems  which  may  be  successfully  used  in  this  way  must  possess 
the  following  characteristics:  (a)  the  active,  transported  chemical  should  have  a 
very  fast  reaction  rate  at  the  active  electrode  surface;  (b)  the  active  chemical 
should  be  reasonably  stable  so  that  no  not  ceable  changrs  in  concentration  occur 
while  the  experiment  proceeds;  (c)  a  supporting  electrolyte  should  be  present  in  a 
concentration  sufficient  to  eliminate  electrostatic  field  effects  in  the  bulk  fluid, 
and  such  that  the  polarization  resistance  at  the  i  surface  is  always  the  controlling 
factor;  (d)  the  rate  of  any  other  electron  trans*.T  reactions  taking  place  must  be 
sufficiently  slow  so  as  not  to  affect  the  primary  reaction  rates  studied. 

Reliable  quality  gauges  of  an  electrochemical  test  system  are  supplied  by  the 
following:  (a)  the  ''  residual  current”,  measured  with  the  experimental  electrochemical 
apparatus  under  test  conditions,  however  with  the  active  species  in  the  solution 
missing,  is  a  measure  of  the  rate  of  extraneous,  unwanted  reactions  taking  place 
concurrently  with  the  rate  which  is  to  be  monitored.  This  residual  current  may  be 
reduced  to  very  low  values  through  the  use  of  clean  soKont  and  chemicals,  and  the 
protection  of  the  solution  from  the  absorption  of  atmospheric  oxygen  and  from 
incident  light  when  necessary;  (b)  the  functional  dependence  of  the  current  on  the 
E.  M.  F.  applied  will  enable  a  conclusion  to  be  reached  as  to  the  suitability  of  the 
system  and  the  activity  of  the  exchange  surface.  When  the  current  measured  (with 
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the  electrode  at  rest,  or  at  any  one,  steady  speed)  is  plotted  versus  the  voltage 
the  limiting  current  should  indeed  be  constant  ov'er  a  substantial  range  cf  E.  M.  F. 

(of  the  order  of  one  volt).  The  level  of  this  current  "plateau"  should  be  unaffected 
by  cycling.  A  deterioration  in  reproducibilit*/  of  the  polarograph  and  a  "sluggishness" 
in  attaining  a  stationary  value  of  the  current  when  convective  flow  or  v'oltage  are 
changed  may  be  an  indication  of  surface  tnsuitability  or  poisoning.  Most  surfaces 
appear  to  deteriorate  after  a  period  of  operation:  this  phenomenon  does  not  appear 
to  be  well  understood. 

The  first  system  to  be  chosen  for  study  was  the  reduction  of  oxygen  on  a 
metallic  gold  surface  in  a  (nominally)  neutral  solution,  with  KCl  (or  KCIO^)  as 
a  supporting  electrolyte,  '^his  system  was  reputed  to  have  worked  satisfactorily 
before  (Dr.  P.  Stonehart);  nevertheless  no  proper  limiting  current  plateau  could 
be  obtained  with  the  test  system  here  'iescribed,  and  reproducibility  of  data  was 
generally  poor.  It  has  the  added  disac:  aniage  that,  with  the  method  of  oxygen 
saturation  recommended  for  this  case  by  the  electrochemists  (simply  bubbling  air 
through  the  solvent)  no  reliable  value  for  o.xygen  concentration  could  ever  be  ob¬ 
tained,  while  at  the  same  time  the  solution  must  have  contained  appreciable  amounts 
of  contaminant  gases. 
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Later  a  more  classical  redox  system.  Fe(CX)  Fe(CN),  was  used. 

t)  D 

This  system  has  been  well  documented  in  tjic  rele\ani  literature  and  property  values 
are  well  known.  Tests  were  performed  with  several  concentrations  of  active 
ions  ,  and  with  different  supporting  electrolytes:  for  all  test  series  up  to  and 
including  XI  both  the  trivalent  and  the  tetravalcnt  salts  in  1/100  normal  concentra¬ 
tion  each  were  used  (Vielstich,  Stonehart),  with  KCl  excess:  for  test  series 

XII  and  XIII  the  concentr adon  of  the  hexacyanoferrates  was  increased  to  5/100 
equivalent  per  litre,  to  reduce  the  importance  of  any  stray  reactions  and  residual 
currents  (Dr.  N.  Ibl,  [23]  ),  in  excess  NaOH.  lastly,  in  test  series  XIV  only  the 
ferri  salt  was  used,  as  all  tests  were  performed  with  the  exchange  surface  electro¬ 
negative.  This  last  composition  of  the  solution  seemed  to  lead  to  the  best  repro¬ 
ducibility  of  results. 

3.  3  Experimental  rig 

3.3-1  General 


Figure  (1)  shows  the  expe r irr  ^ntal  equipment  schematically;  figure  (2)  is  an 


SCHEMATIC  DRAWING  OF  TEST  APPARATU 


FIGURE  2,  DIAGRAM  OF  ELECTRICAL  CONNECTIONS. 
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electrical  circuit  diagram  while  the  plate,  figure  (3),  shows  the  assembled  equip¬ 
ment. 

The  apparatus  .isists  essentially  of  a  cylindrical  vessel  of  soda  lime  glass 
17-1/2  inches  in  diameter  and  18  inches  deep.  A  "Perspex"  cover  carries  a  drive 
mechanism  to  which  a  rotating  arm  assembly  is  attached.  A  cylindrical  electrode 
carrier  of  1/8  to  1/4  inch  diameter  is  mounted  on  the  rotating  arm,  and  the  latter 
is  driven  in  a  circular  motion  by  the  drive  mechanism.  Thus  the  test  electrode 
moves  along  the  circumference  of  a  circ’e  of  about  4.  62  inches  radius.  The  glass 
container  is  filled  with  electrolytic  Huid  to  a  height  of  about  5-3/4  ir.''hes. 

A  constant  voltage  transformer  supplies  an  induction  motor  which  m  turn  drives 
the  rotating  arm  assembly  through  assorted  step-down  devices. 

The  D.  C.  to  the  test  bath  is  supplied  from  a  potentiostat:  the  test  surface  is 
usually  acting  as  the  cathode.  The  anode  is  a  platinum  foil  "counterelectrode". 
Current  '»  taken  off  through  a  trailing  wire  connected  to  the  test  electrode.  This 
method  of  current  take*cff  was  chosen  for  its  simplicity  over  the  more  elegant  and 
orthod(5x  mercury  pool  slipring  or  sintered  carbon  brushes:  thus  also  the  number 
of  items  requiring  developn  ent  was  reduced  by  one. 

Much,  time  was  spent  on  insuring  that  both  the  test  vessel  and  its  cover  were 
mounted  horizontally.  To  reduce  free  surface  effects  the  active  test  e’  ctrode 
area  extended  only  a  comparatively  short  vertical  distance, deep  inside  the  body 
of  the  fluid.  Finally,  to  insure  a  steady  rate  of  rotation  of  the  electrode  carrier 
the  driving  motor  used  was  of  ample  size,  and  the  final  step  down  pinion  was 
mounted  on  ball  bearings.  Anti- vibration  mounting  was  used  throughout. 

During  initial  tests  at  very  low  transfer  rates  and  speeds  it  was  found  that 
strav  vibration  and  cantilever  bending  resonance  of  the  test  electrode  carrier 
were  on  occasion  increasing  transfer  rates  notably.  It  took  mv'-h  effort  to 
cT-minate  these  secondary  effects,  but  the  phenomenon  pointed  towards  the 
potential  usefulness  of  the  experimental  method  for  the  measurement  of  insta- 
tionary  mass  transfer  processes.  In  this  respect  the  present  technique  has  the 
advantage  over  heat  transfer  measurements  in  that,  when  used  with  a  fluid  of 
high  Schmidt  number  (which  it  almost  invariably  is),  there  are  almost  no  time 
lags  of  transfer  rate  associated  with  a  ''hange  in  convective  rate.  Thus  the 
resolving  power  of  this  techni  lue  should  be  potentially  higher  than  that  of  heat 
transfer  experiments. 
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Froude  number 


13 


FIGURE  3  experimental  EQUIPMENT 
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3,3-2  Electrodeg 

Initially  a  solid,  metal  cylindrical  electrode  {1/8"  diameter)  was  used.  The 
metal  cylinder  was  coated  with  an  insulating  lacquer  and  a  slit  about  1  cm  long  by 
0.  15  cm  wide  was  cut  in  the  lacquer.  This  then  formed  an  active  electrode  area 

the  size  of  which  was  known  with  reasonable  accuracy.  However,  once  the  electrode 

was 

carrier  /  assembled  only  a  single  angular  orientation  could  be  tested  at  a  time.  The 
test  surfaces  were  usually  made  of  one  of  the  noble  metals,  so  that  this  mode  of 
construction  was  also  rather  expensive.  Finally,  none  of  the  insulating  lacquers 
used  would  last  very  long  under  actual  test  conditions. 

Later  "Bakelite"  and  "Ebonitv"  carriers  of  various  diameters  were  developed. 
These  had  axial  grov^es  milled,  into  which  wires  of  the  test  electrode  material  were 
glued  with  epoxy  resin.  An  actsv'e  electrode  area  was  then  obtained  by  grinding 
a  length  of  about  3/4  inches  of  <  arricr  length  in  such  a  manner  that  a  fresh  metal 
surface  lay  bare  on  the  cylinder  envelope.  These  test  pieces  were  inexpensive  and 
sturdy,  and  the  number  of  angular  orientations  which  could  be  tested  wiihout  dis- 
nn.antling  the  rotating  arm  assembly  depended  only  on  the  number  of  axial  groves 
which  could  be  conveniently  milled  in  the  carrier.  There  was  however  one  disad¬ 
vantage  in  that  the  active  test  area  was  not  known  accurately  and  thus  calibration 
tests  had  to  be  run. 

.Electrode  materials  used  were  annealed  gold,  platinu.m  and  nickel  wires.  The 
wires  were  of  0.010"  to  0.025"  diameter.  The  counterelectrode  was  an  annealed 
platinum  foil  on  all  test  runs. 

One  of  the  objects  of  the  study  here  reported  was  the  study  of  flow  conditions 
in  regions  of  separated  How  behind  bluff  bodies  immersed  in  a  flow.  In  particular, 
it  was  hoped  that  the  degree  of  stabilization  achieveable  by  splitter  plates  placed 
downstream  of  the  electrode  carrier  could  be  demonstrated.  Accordingly  several 
carriers  with  attached  splitters  were  produced,  at  first  with  stainless  steel  vanes 
coated  with  an  epoxy  resin,  later  with  platinum  vanes.  Electrode  carrier  cot;struc- 
tion  and  a  carrier  with  vane  are  shown  in  the  plates,  figures  (4)  and  (5). 
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3.3*3  Preparation  of  the  test  electrode  surface 

For  the  first  test  series^using  the  oxygen  reduction  reaction  on  a  gold  svirface^ 
no  preparation  apart  from  cleaning  of  the  gold  surface  with  an  organic  solvent 
was  undertaken  {Dr.  P.  Stonehart).  For  the  last  of  the  hexacyanofer rate  series 
the  activation  procedure  used  by  the  Zlirich  tean  (Ibl,  [^23"]),  was  employed; 
after  each  test  series  the  surface  of  the  electrode  was  first  ground  to  a  smooth 
finish,  whereupon  electrode  and  carrier  were  cleaned  briefly  in  trichloroethylene 
and  rinsed  in  distilled  water.  The  test  electrode  was  then  electrically  cycled  for 
about  5  minutes  in  a  solution  of  about  10  ppw  KCN,  30  ppw  Na^CO^,  30  ppw  KOH 
in  500  ppw  of  aqua.  The  test  electrotle  was  connected  cathodically  and  the  applied 
voltage  was  raised  until  gas  evolution  took  place.  The  electrode  thus  prepared 
was  rinsed  and  immersed  in  the  test  bath. 

The  counterelectrode  i  nd  the  free  convection  foils  were  similarily  prepared. 

3. 4  Measuring  equipmei.. 

3.4- 1  The  referer.ee  electrode 

A  standard  mercury  pool  type  reference  electrode  was  used  in  order  to  measure 
the  potential  difference  between  the  test  electrode  surface  and  the  solution.  For 
the  first  test  series  a  calomel  solution  in  KCl  was  used  for  contact  with  the 
mercury,  later  HgO  (mercuric  oxide)  in  NaOH.  In  order  tc  reduce  the  diffusion 
of  mercury  salts  into  the  test  bath  a  special  type  of  test  electrode  holder  was 
developed,  figure  (6). 

The  voltage  between  reference  electrode  and  test  piece  was  measured  at  no 
current,  by  an  instrument  with  a  100  Mfl  input  impedance. 

3.4- 2  The  current  was  measured  by  inserting  a  low  impedance  galvanometer 
in  the  return  from  the  test  electrode.  The  period  of  rotation  was  measured  by 
electrical  timers  which  were  first  calibrated. 

The  accuracy  of  the  instrumentation  used  finally  exceeded  by  far  the  accuracy 
and  reproducibility  of  the  system  as  a  whole,  so  that  a  detailed  analysis  of 
measurement  error  did  not  seem  to  be  warranted. 

3.  5  Data  reduction 


3.5-1  Theory  of  electrochemical  process 
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FIGURE  6  REFERENCE  FLECTRODE  AND  ELECTRODE  HOLDER 
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Assuming  a  single  reaction  at  the  test  electrode  k;urface,  the  transfer- rate  is 
given  by  the  following  expression; 


TR  = 


F  X  N. 


..{3.  1) 


where  i”  is  the  interface  current  density;  F  is  the  Faraday  constant  and  is 


the  reaction  valence.  The  mass-transfer  coefficient  is  obtained  as  follows: 


H  = 


c  -  c 
00  o 


..(3.2) 


where  Cq^  is  the  concentration  of  the  active  species  in  the  bulk  fluid  and  co  is 
the  interface  concentration.  This  latter  value  is  commonly  assumed  to  be  zero 
for  reactions  proceeding  at  a  fast  rate.  An  expression  for  the  mean  Sherwood 
number  is  obtained  by  inserting  (3.  1)  into  (3.2): 


Sh  = 


A  •  (Cqo  -Cq  )NvF^ 


...(3.3) 


where  I  is  the  current  measured;  I  is  a  reference  length;  A  is  interfacial  area 
and  %  is  the  diffusion  constant. 

3.  5-2  Calibration  tests 

For  most  of  the  test  conditions  neither  the  interfacial  area  nor  the  diffusion 
constant  were  known  with  great  accuracy  a  priori.  Independent  calibration  tests 
were  therefore  run  as  follows.  The  diffusion  constant  was  determined  by  meas¬ 
uring  the  natural  convection  mass  transfer  rates  on  platinum  foils  of  accurately 
known  dimensions,  suspended  vertically  in  the  fluid.  For  this  case  the  mean 
dimensionless  transfer  coefficient  at  large  Sc  number  is  given  by, 


Sh  =  0.  653  (Gr-Sc)^ 

where  the  reference  length  I  is  plate  height. 
Inserting  (3.  3)  for  Sh  and  solving  for  ^  , 


^  =  0-7009  ^  ) 


(3.4) 


(3.  5) 


where  b  is  plate  width.  Or,  in  dimensionless  forms. 
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V 


-Vj 

=  0.  7009  Gr  ’ 


b  F  (c^-c^) 


(3.  b) 


U5U?Uy,  at  least  two  tests  were  run,  using  plates  with  different  dimensions. 
Values  of  Gr  obtaining  in  practice  were  in  the  range  4  x  10  to  3  x  10^  In 
order  for  the  asymptotic,  boundary  layer  type  theory  outlined  above  to  be  appl: 
cable,  we  must  have. 


Sc  >  10' 


5  X  10  <  (Gr  •  Sc)  <  10 


The  upper  limit  on  Gr'Sc  is  introduced  by  the  onset  of  transition  to  a  turbulent 
boundary  layer.  Also,  in  order  to  reduce  three-dimensional  flow  effects,  the 
ratio:  foil  height  to  foil  width,  rt  Ih,  should  be  smaller  than  I  M. 

In  order  to  be  able  to  evaluate  the  Grasshof  number  at  test  conditions,  a 
knowledge  of  the  change  in  density  taking  place  at  the  test  electrode  surf'’co  i 
necessary.  The  relevant  data  for  the  system  here  used,  but  at  different  conct  n- 
tra'^ions,  have  fortunately  already  been  previously  reported  ^24']  .  The  accuracy 
of  which  the  experiment  described  here  was  capable  did  not  warrant  an  indepondt  nt 
detailed  check  of  these  data,  which  can  be  done  by  using  tiie  method  of  Fouad 
and  Ibl  C25  J ,  and  interpolation  data  between  the  reported  test  points  C-^^wert  used. 

The  kinematic  viscosity  of  solutions  of  the  composition  and  concentrations 
used  in  the  present  series  of  experiments  has  prev  iously  been  reported  (c.f. 

C20I,  p.  530),  but  these  were  in  part  repeated  independently.  ' 

Chce  the  Schmidt  number  v/as  known,  the  test  electrode  area  was  obtained  by 
assuming  equation  (2.  b)  to  be  applicable  at  the  front  stagnation  point  when  the 
test  electrode  is  measured  in  forced  coneection.  Insertion  of  - 


fi  <x) 


3.48b5^’  =  3.4865X 


(3.  7) 


x«  =  0 


into  (2.  6),  (2.  6a)  and  comparison  with  eq.  (3.  3)  yields, 

4 


(x)  =  g  /3 


A  =  1.069 


I.  R 


C  F£) 


Re  Sc 


(3.  8) 


Values  found r)/-  1 .  32 5  cent! stokes  at  23  C;  1.267  at  25;  1.2  12  at  27;^  -  1.090  at  23,9  C 
Grnssmann  el  al.  24  reported  a  value  of  1^*1 .  2 1  centistokc s .  See  also  Appendix  III. 
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3.  5-3  The  direct  measurement  of  wall  shear  stress 

Once  both  the  diffusion  coefficient  and  the  interfacial  area  hav'e  been  determined 
the  local  wall  shear  rate  at  locations  away  from  the  front  stagnation  point  may  be 
obtained  directly  by  eliniinating  Sh  from  equations  (2.  6)  and  (3.3)  for  the  constant 
wall  concentration  case: 


T 


_A^ 
3  Sc 


.  .  .(3.9) 


3,.  f)  Results 

The  graph  figure  (7)  shows  an  example  of  a  polarograph  obtained  during  one 
of  the  first  test  series  with  the  oxygen  on  gold  reaction.  It  will  readily  be 
realized  that  there  is  no  real  limiting  current  "plateau"  and  that  the  effect  of 
cycling  hysteresis  is  very  notable.  This  latter  phenomenon  may  hav'e  been  due 
both  to  continuous  change  in  oxygen  saturation,  or  rather  supersaturation,  on 
the  one  hand  a.nd  surface  activation,  later  counteracted  by  surface  poisoning,  on 
the  other. 

The  graph,  figure  (8)  shows  an  example  of  polarograph  obtainable  with 
hexacyanofer rates  redox  reactions  and  with  a  reasonably  clean  technique.  When 
quite  fresh,  electrode  activ'e  areas  may  yield  limiting  current  "plateaux"  which 
extend  virtually  from  +0.25  to  -1.0  volts  without  a  noticeable  change  in  the 
value  of  the  current  passed. 

The  graph,  figure  (of  gives  transfer  rates  for  a  system  of  which  graph  (8) 
represented  the  polarograph.  It  is  seen  that  the  front  stagnation  point  transfer¬ 
rate  values  (marked  FSP  on  the  graph)  follow  the  theoretically  expected  line 
within  a  reasonable  amount  of  scatter.  Surprisingly,  though,  the  values  for  tho 
rear  stagnation  point  (marked  RSP)  also  fall  on  a  line  with  slope  1/^  n  the  log- 
log  plot.  Values  for  an  angle  of  tt/2  fall  above  those  for  0°  and  for  tt,  as  expected 
from  the  theory  (equation  2.6  above). 

In  the  graph  figure  (10)  results  are  giv'en  for  the  case  of  cylindrical  electrode 


’^'Figure  (9)  is  not  directly  comparable  with  the  following  figures  due  to  the 
different  concentrations  of  chemicals  used. 
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carriers  with  and  without  a  stabilizing  vane  attached  to  their  downstream  side. 

The  plate  figure  (5)  showed  an  electrode  carrier  ith  attached  stabilizing  vane. 

It  is  seen  from  graph  (10)  that  the  level  of  transfer  densities  is  reduced  when 
the  electrode  carrier  with  a  vane  s  tested,  as  compared  to  the  same  carrier 
when  the  vane  was  cut  off.  The  character  of  the  curves  is,  however,  not  altered. 

Figure  (11)  illustrates  the  results  when  the  electrode  carrier  w’ith  vane 
attached  is  rotated  in  the  opposite  direction,  so  that  the  \  ane  faces  upstrearrj. 
Predictably, transfer  rates  are  reduced, at  what  is  now  the  F.  S.  P.,  in  the  low 
Reynolds  number  range,  as  the  F.  S,  P.  electrode  is  now  bathed  in  the  viscous 
wake  from  the  upstream  splitter  vane  and  local  wall  shear- rate  values  are 
reduced.  At  high  Re,  however,  transpor  rates  are  incieased:  this  may  be 
attributed  to  the  fact  that  the  upstream  splitter  vane  sustains  a  non-zero  a.-gh 
of  incidence  to  the  flow,  and  thus  a  sizeable  \-elocity  from  the  pressure  side 
of  the  pla  ,e  to  the  suction  side  may  sweep  past  the  electrode  surface. 

The  graph  figure  (12)  shows  complete  results  of  test  series  XIV  giving 
v'alues  of  transfer  coefficients  at  various  Reynolds  numbers,  with  the  angle 
substended  to  the  flow  as  a  parameter,  The  small  amount  of  scatter  obtained 
is  notable. 

3.6-1  Comparison  of  the  experimental  results  with  the  theory  of  para.  2.  !-l 

In  the  graph,  figure  (13)  the  results  given  in  the  previous  figures  are  re¬ 
drawn  so  as  to  give  the  ratio  of  Sherwood  number  to  (Sherwood  number  with 
the  mass  transfer  strip  at  the  front  stagnation  line),  plotted  versus  the  circum¬ 
ferential  angle  substended  by  the  median  of  the  strip,  and  with  the  Reynolds 
number  as  a  parameter.  A’so  drawn  in  the  graph  is  the  result  obtained  from 
equation  (2.  5)  based  on  the  follov  ing: 

An  experimentally  determined  velocity  distribution  aroitnd  a  cylinder  in 
stabilized  laminar  channel  flow  has  been  obtained  by  Grove  [29]  .  and  hi.'s  data 
may  be  correlated  by  the  following,  Acrivos  C3o  J  , 

u^  =  3.  758^-^-  2.84^”*  +  0.674^''  ...(3.10) 

Here  the  velocity  u  has  been  rendered  dimensionless  by  using  as  reference 
velocity  the  value  obtained  %/ithin  the  channel  at  the  location  of  the  cylinder 
before  the  latter's  insertion.  In  view  of  the  fact  that  the  blockage  ratio  of  the 
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channel  was  appreciable,  ( 3.  10)  cannot,  of  course,  be  t' ken  to  represent  a  good 
correlation  for  the  widely  differing  case  here  considered.  Another  esthetical 
flaw  is  that  (3.  10)  is  not  symmetrical  with  respect  to^.  It  was  used  neverthe¬ 
less,  as  there  are  few  reliable  experimental  data  available. 

From  (3.  10)  the  local  shear  stress  at  the  wall  ma  '  be  calculated  using  one 
of  the  approximate  methods.  We  used  here  that  of  Acrivos  ri3jand  obtained. 


/3(p  =  I  [4.97  -  4.84yl'+  1,157^''1  ...n.ll) 

The  v'alues  of  thus  obtained  may  be  introv^uced  into  equation  (2.  5)  together 
with  the  requisite  source  -  strip  width  (11.  b  degrees  in  the  present  case)  to  give 
an  estimated  curve  of  the  circumferential  distribution  of  mass  transfer  coef¬ 
ficients.  It  is  seen  from  the  graph,  figure  (13),  that  there  is  fair  agreement 
over  that  part  of  the  cylinder  circumference  where  such  an  empirical  relation 
could  be  expected  to  yield  reasonable  results,  id  est  not  too  close  to  the  boundary 
layer  separation  point.  The  results  plotted  in  graph  (13)  thus  gi\'e  us  some 
confidence  in  using  the  experimental!/  determined  mass  transfer  data  in  order 
to  estimate  wall  shear  rates  in  the  range  ^2*^  to  180^^  where  no  other  direct 
measurements  have  been  available  heretofore. 

3.6-2  Critique  of  method  and  suggested  improvements 

The  main  difficulties  experienced  were  electrochemical  in  nature:  (a)  the 
diffusion  coefficient  found  from  our  free  cori'.ection  experiments  is  rather  low- 
compared  o  that  reported  by  other  investigators.  The  v'alue  found  here  was 
3.  1  X  10  ^  cm^  sec^  in  test  ser-es  XIII;  3.  58  x  10  ^  at  21*^0  and  4.  13  x  10 
at  25.  3°C,  in  test  series  XIV.  The  value  reported  by  Grassmann  and  coworkers 
is  5.8  X  lO'^^  at  25°df^b)  after  a  very  short  time  of  operation  the  response  time 
of  the  system  to  altered  conditions  of  speed  or  driving  force  (E.  M.  F.  )  became 
very  large  indeed.  Once  the  surface  had  thus  deteriorated,  values  of  current 
measured  v  nder  any  given  conditions  were  invariably  higher  than  the  true 
steady-state  values  obtained  with  a  fresh  surface;  (c)  it  would  appear  likely  that 
in  order  to  obtain  consistently  reproducible  experimental  data  a  continuous 
activation  of  the  test  electrode  surface  by  sup  -position  of  an  A.  C.  component 
on  the  D.  C.  pa  ;ed  through  the  system  is  necessary. 

(*)  ~ 

The  use  ol  i  iv  absolute  -.-alue  for  the  second  term  v.'ould  render  the  function 

non  analytical.  Si  milar  rema;  ks  apply  to  equation  (3.  11). 
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tially  two-dimensional  flow 

As  pointed  oui  in  section  3.3-1  the  electrochemical  system  described 
measures  transfer  -  rate  oscillations,  with  amplitvidc  and  frequency  distortions 
v/hich  can  be  made  arbitrarily  small  by  choosing  suitable  measurement  instru¬ 
mentation.  Thus  it  is  singularily  well  suited  to  the  measurement  of  the  vortex 
shedding  frequency  behind  blunt  objects  placed  in  a  stream,  when  the  monitor¬ 
ing  electrode  is  placed  near  the  rear  stagnation  point:  the  oscillatory  frequency 
in  the  current  measured  corresponds  to  a  synchronous  oscillation  in  the  local 
shear-rate  at  the  wall. 

It  is  shown  in  classical  work  that  the  oscillatory  frequency  can  be  correlated 
with  the  linear  dimensions  of  the  obstacle,  the  fluid  \-elocity  and  fluid  kinematic 
viscosity  in  a  unique  way: 

S  =  F  (Re)  ...  (3.  12) 

where  S  is  the  Srouhal  number,  S  =  f.  1/U  ,  with  f  equal  to  the  frequency; 

(ZD 

a  characteristic  length  (usually  the  greatest  extent  of  the  body  in  a  direction 
perpendicular  to  that  of  the  flow  velocity  at  great  distance  from  the  body, 
and  Re  is  the  Reynolds  number.  The  functional  relationship  (3.9)  cannot  be 
determined  from  analysis  and  must  be  found  from  tests.  However  F  is  very 
little  dependent  on  the  detailed  cross-sectional  shape  of  the  blunt  body  tested. 

In  the  graph,  figure  (15)  measured  values  of  frequency  (obtained  by  timing 
the  oscillation  of  the  current  observed  on  the  galvanometer  in  the  test  electrode 
circuit)  are  compared  to  data  reported  by  Roshko,  [263,  1  >  bhair,  [28^  ■ 

Shair's  data  are  apparently  not  very  reliable  >[30  '  as  in  his  tests  the  v'elocities 
were  incorrectly  measured:  however,  in  conjunction  with  Roshko's  very 
careful  measurements  and  the  data  reported  here,  they  seem  to  indicate  that  at 
any  given  Reynolds  number  at  which  instabiUty  does  indeed  take  place,  the  Strouhal 
number  is  the  higher  the  larger  the  channel  blockage  ratio.  This  result  cannot 
be  taken  too  seriously  without  much  additional  careful  testing. 

Another  interesting  point  apparent  from  graph  (15)  is  that  while  the  stabilizing 
vane  did  not  suppress  instability  of  the  wake  within  the  range  of  Reynolds  numbers 
where  it  was  expected  to  dc  so,  it  nevertheless  reduced  gross  turbulence  somewhat 
and  thus  gave  rise  to  lower  Strouhal  numbers.  The  same  effect  was  found  by  Roshko, 


0.30 


ncr eased  25%  to  allow  for  incorrect  doto.resuNson  this  graph  show  sliding  scole should  be  used. 
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4 .  Experimental:  Heat  Transfer  Studies 
4.  1  Introduction 

It  was  pointed  out  in  sections  3.  1  through  3.  2  above,  that  the  initial 
expectations  were  that  the  electrochemical  mass  transfer  studies  vvoulH 
yield  rapid  and  accurate  results  with  little  expenditure  and  complica¬ 
tion  From  the  expose  in  section  3.6  it  will  howecer  have  become  clear 
that  the  accuracy  and  the  r eproductibility  of  these  experiments  were  rather 
■vorse  than  one  was  led  by  elect rochemists  to  believe.  Also,  no  conclusive 
evidence  of  the  degree  of  stabilization  achiecable  with  the  splitter  plate, 
and  the  influence  of  this  stabilization  on  RSP  mass  transfer  rates  for  flows 
with  very  low  blockage  ratios  was  forthcoming.  So  it  was  decided  to  run 
some  equivalent  heat  transfer  studies. 

A  closed  circulation  laminar-flow  oil  tunnel  was  available  in  the 
Department  and  has  alreadv  been  described  previously  [l4.,  [is],  Cl?], 
[28],  C29J.  l31J.  C32*.  It  IS  apparent  from  those  previous  studies  that 

this  oil  tunnel  may  not  be  an  ideal  research  tool  for  the  study  here  proposed 
due  to  some  uncertainty  in  the  interpretation  of  previous  experimental  data 
(compare  C23J  through  [32J).  Enough  is  nevertheless  known 

about  the  variation  of  the  rates  of  shear  on  the  wall  of  a  test  cylinder  placed 
in  that  tunnel,  at  least  up  to  the  point  of  boundary  layer  separation.  Thus 
data  to  be  measured  in  that  region  may  be  used  for  calibration  purposes. 


4.  2  Construction  of  heat  transfer  cylinder 

The  diagram  figure  (16)  and  plate  (17)  show  the  construction  of  the  heat 

transfer  cylinder  with  a  strip  heat  source.  Outer  diameter  was  3/4  inches 

nominal.  It  will  be  seen  that  this  cylinder  is  of  simple,  rugged  and  cheap 

(«) 

construction.  The  heat.ng  surface  consisted  ot  a  carefully  rolled  copper 
strip,  1/16  inch  wide  and  0.  009  inches  thick.  It  was  embedd'^d  in  a  grove 
cut  parallel  to  t*  •  cylinder  generating  line.  The  copper  h<->atlng  strip  w'as 

F) 


The  help  of  tr.e  Mechanical  Workshop  of  the  Chemistry  Department  in 
the  const nu  tion  of  the  heat  transfer  cylinder  is  gratefully  acknowledged. 
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FIGURE  17  TEST  HEAT  TRANSFER  CYLiMDFR 
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FIGURE  15  /^55EMgLeP  TEST  CYLINDER 
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"short-circuited"  over  the  terminal  1  inch  near  the  channel  confining 
walls,  30  that  appreciable  heating  took  place  only  over  the  central  6 
inches  of  the  channel.  This  was  done  in  order  to  avoid  end  effects 
through  the  heating  of  the  rather  thick  wall  boundary  layers.  The  short* 
circuiting  was  carried  out  by  having  the  h  ating  strip  at  5  times  its  "cen- 
tral'lhi cknes s ,  and  d  epening  the  grove  holding  the  strip  corresponding¬ 
ly. 

Three  thermocouples  were  welded  in  the  standard  manner  using 
copper  constantan  32  BWG,  fiberglass  arid  enamel  insulated  wires.  The 
cleaned  junction  w-as  then  silver  soldered  to  the  underside  of  the  strip 
at  locations  dividing  the  heated  length  into  quarters.  The  thermocouple 
wires  were  led  to  the  central  bore  of  the  cylinder  and  thence  to  the 
measurement  apparatus.  Finally  the  heating  strip  was  glued  into  place 
using  an  epoxy  resin  ("Araldite  '). 

Three  other  thermocouples  w’ere  prepared  in  the  same  manner  and 
glued  into  countersunk  holes  at  various  circumfeiential  locations  on  the 
cylinder  er.velope.  These  latter  thermocouples  served  for  the  determin¬ 
ation  of  the  free-stream  temperature. 

At  each  end  of  the  cylinder  a  heavy  brass  ring  was  shrunk  with  a 
slight  interference  fit  over  the  otherwise  assembled  cylinder,  and  the 
heating  strip  was  silver  soldered  to  the  rings.  Thus  a  good  electrical 
contact  was  assured.  The  brass  ring  serv  ed  both  for  supply  of  the 
rather  heavy  low  voltage  current  and  for  carryine  "0"  rings  assuring  a 
close  fit  of  the  cylinder  in  the  channel  ’walls. 


4.  3  Measurement  Technique 

Current  was  supplied  fre-m  a  high  current  AC  source  '.  The  vol¬ 
tage  drop  over  the  heat.ng  strip  was  measured  by  coimecting  a  ver’/  high 
impedance  AC  voltmeter  to  the  brass  current  supply  rings  described 

{*)  Thanks  are  due  to  the  Department  of  Mechanical  Engineei  ing  at  Stan¬ 
ford  for  the  loan  of  the  current  source. 
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above.  On  a  later  build  separate  leads  were  soldered  to  the  heating  strip, 
so  that  only  the  voltage  drop  over  the  central  6  inches  of  strip  was  mea¬ 
sured,  making  for  greater  measurement  precision. 

The  resistance  of  the  heating  strip  wai  deter*  ined  by  separate  calibration. 
A  constant  DC  current  delivered  by  a  potentiostat  was  measured  by  n  pre¬ 
cision  H.  P.  milliammeter  and  passed  through  the  strip.  The  voltage  drop 
was  measured  with  the  aid  of  an  electronic  microvoltmeter  (COHU- 
KAYLAB  model  203).  Values  found;  7.41  milliohms  for  the  active  6  inches 
of  heating  strip,  and  12.9  milliohms  between  the  supply  clips,  at  an  ambient 
temperature  of  22.1  °C. 

Temperatures  were  measured  in  the  usual  way,  using  a  separate  cold 
junction  connexion  for  each  thermocouple.  A  Leeds  and  Northrupp  K3  po¬ 
tentiometer  was  used,  which  enabled  readings  to  be  taken  to  1  micrcv'clt . 

The  tunnel  test  variables  were  measured  as  described  in  28  and  29  . 
and  will  be  described  in  detail  in 


4 .  4  Result ^ 

In  contradistinction  to  the  mass  transfer  experiments  described  in  chapter 
(3)  above,  the  heat  transfer  experiments  proceeded  without  fuss  or  running -ir. 
difficulties.  Indeed,  the  very  first  calibration  rvin  yielded  useful,  consistent 
and  reproducible-  results,  figures  (19)  and  (20).  The  test  apparatus  proved 
quite  sensitive,  and  such  detailed  features  as  the  slight  variation  of  transfer 
coefficients  inside  the  stabilized  wake  bubble  with  change  in  attitude  angle  or 
in  the  positior  of  the  splitter  vane  could  be  clearly  discerned.  It  is  seer,  from 
the  graph,  figure  (20),  that  the  agreement  with  the  theory  developed,  using 
Grove's  empirical  velocity  correlation,  is  surprisingly  good.  A  somewhat  un 

(*) 

I  am  indebted  to  Mr.  l>.  Snowden  for  his  cooperation  in  the  use  of  the  oil 
tunnel  and  lor  his  help  in  running  the  equ  nment. 
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expected  result  is  the  notable  variation  of  transfer  coefficients  as  the 
the  channel  is  traversed  crosswise.  This  would  point  towards  notable 
wall  effect?.  The  central  thermocouple  does  indeed  read  temperatures 
which  are  considerably  higher  than  those  read  by  the  two  lateral  ones. 
Doubtlessly,  thermal  conduction  longitudinally  along  the  strip  could 
also  have  been  of  some  importance,  though  it  is  believed  that  this 
latter  effect  cannot  have  been  v'ery  important. 

The  cylinder  with  the  strip  type  heat  transfer  surface  was  turned 
over  to  Mr.  Snowden  for  further  teats,  which  will  be  reported  upon 

in  [32] . 
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Appendix  I. 

Exparnion  of  u  (x,  y)  near  a  flow  boundary  at  which  a  no- slip  condition  applies. 
At  any  solid  boundary  of  a  flow  the  following  conditions  apply: 


u 


1  =0 


.(I-l) 


where  and  are  the  dimensional  velocities  in  the  direction  (parallel 
to  the  boundary)  and  the  y^  direction  (perpendicular  to  the  boundary)  res¬ 
pectively.  The  momentum  equation  at  the  boundary  thus  degenerates  to, 


9^u^ 


8x^  y*  =0  ■  ° 

Assuming  that  u^  is  an  analytic  function  of  the  ordinate  y.^^ 


.  .  .(1-2' 


Uj.(x^,y^) 


X  =  constant 


=(S  <>'1 ) 


then, 


Now,  from  (I-  1) 


and  from  (1-2), 


=  /(O)  +  +/"(0)  fr  +  •• 

^  (0)  =  ;  also. 

/•'  (0)  =  9“  '  To(’‘'/yW 

^■”(0)  =  0 


/T- 


(1-4) 


We  shall  insert  these  results  into  the  momentum  equation  so  as  to  obtain 
the  higher  order  terms.  Integrating  twice  we  find  that  the  expansion  for 
correct  to  the  5th  order  terms  is, 

„  1  -  ‘is.  V  .L  ^  (A  r  ±5  .  Ae  \ 

I  i'Vi  ^  yt  ^  dx  dx^**  24^  *  dx^  dy^-^j 

+  ijL,  (1  +  An ^  + 

bOyLfc  '  2  dx  ^5.  ^  dx^ 


...(1-5) 


If  the  energy  equation  is  to  be  solved  for  fluids  of  very  large  Prandtl  numbers, 
the  expression  for  u^  may  be  linearized,  retaining  only  the  first  term  of  the 
expansion,  as  the  thermal  boundary  layer  is  very  thin  compared  to  the  morr.en- 
tum  boundary  layer.  This  is  essentially  the  Leveque  assumption  Qb,  12^. 
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Curvature  and  thickness  effects  are  neglected.  The  conduction  equation  for 
the  heating  strip  then  reduces  to, 


^  (KJ9)  )  +  q,'"  (0)  -  =  0  ...(II-l) 

where  0  is  the  excess  temperature: 

is  the  conductivity  of  the  solid  wall  material; 

H  is  the  surface  convection  coefficient  and  ^  is  strip  thickness; 

q  '"is  heat  source  strength  per  unit  volume  of  material;  for  electrical 

8 

heating, 

q^'"  =  CV*/(^^L^)  ...(II-2) 


where  V  is  the  applied  electrical  potential  difference  over  a  strip  of  length  1 
and  electrical  resistivity  /«  (9);  C  is  a  dimensional  constant  (equal  to  0.24  if 
q'"  is  in  cgs  units,  V  in  volts  and  ^  in  ohms  cm). 

The  appropriate  boundary  conditions  are, 

‘  X  =  X*  and  tc.  ~  ■¥  ® 


We  shall  now  assume  all  physical  properties  to  be  independent  of  temperature. 
Equation  (I-l)  may  be  written  in  dimensionless  form  thus. 


h 

V 

E 

c 

a 


a 

v; 

u; 

II.  1  F 
tl 
tl 
e 
e 
tl 


T 

Ii 

a 


0 


...(11-4) 
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Here  ^  -•  x/r,  w-th  r  a  reference  length;  is  a  dimensionless  temperature 

=  9/M 

where  M  is  the  source  strength,  M  =  q  r  /K  ,  and  Bi  is  a  Biot  number, 

H  r^  ®  ® 

Bi  =  ~ —  X  -p—  .  Further  introducing  the  asymptotic  value  of  the  heat  transfer 

0 

coefficient  for  large  Pr  numbers,  equation  (2.3);  transforming  the  coordinates 


as  follows, 


x:  = 


T  = 


f 

K  fluid 


Re'*Pr’^  9 


K  solid  ^ 
and  introducing  a  parameter, 


Tc  -  T 

^O  00 


C  _  K'^solid  ^  Tsq-Tqo 
^  ■  K  fluid  rS  q 

s 


9V(I)  _L 

Re'iPr'3 


(II-5) 


we  finally  obtain, 

et 

^  '  XX 


dT,  (X)  _ 


=  0 


{II- 5a) 


(11-6) 


with  boundary  conditions  T'(0)  =  0  and  T'  (I)  =  0- 

such 

1  For  almost  isothermal  heating  of  a  metal  strip/as  in  the  diagram  figure  I-l, 
the  numerical  value  of  the  parameter  £  is  very  large  indeed.  If  we  introduce 
the  additional  simplification  of  replacing,  over  the  strip  width  in  the  flow  dir¬ 
ection,  yQ  QQ  by  its  average  value,  then  a  regular  perturbation  solution  of 
equation  (1-6)  is  readily  obtained.  This  yields  an  expression  for  the  tempera¬ 
ture  distribution  at  any  location  in  ascending  powers  of  ^  : 


+  0 


Ir 


T.  .J  (I  -^T^) 


(II-7) 


The  constant  has  to  be  determined  from  a  heat  balance  calculation. 

In  the  practical  case  of  the  heat  transfer  cylinder  with  a  copper  strip  described 
above  the  numerical  value  of  the  parameter  £  was  about  21. 
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n.2 

We  ahall  now  consider  the  case  of  a  metal  strip  which  forms  an  almost 
constant  flux  surface:  i.e.  the  parameter  £  is  now  very  small.  While  in  the 
previous  case  a  passage  to  the  limit  £->oo  did  not  decrease  the  order  of  the 
differential  equation  im'-olved,  and  hence  a  regular  perturbation  scheme  could 
be  evolved,  in  the  case  to  be  now  analyzed  the  limit  £-»0  would  decrease  the 
order  of  the  equation:  we  shall  therefore  expect  to  have  to  evolve  a  singular 
perturbation  scheme.  This  is,  however,  not  straightforward  in  the  present 
case. 


Considering  the  actual  physical  configuration,  the  importance  of  the  second 
derivative  in  the  equation  (the  conduction  term)  is  significant  only  near  the 
boundaries  of  the  strip,  where  this  tei  m  will  have  to  insure  the  fulfillment  of 
the  boundary  conditions.  Over  the  rest  of  the  strip  a  "core  solution"  will  be 
operative  and  conduction  v/ill  be  unimportant. 

Proceeding  as  in  the  case  supra,  our  basic  equation  is  still  (II-4),  however 
now  with  a  very  large  value  of  Bi.  Replacing  H  from  equation  (2.  12)  we  now 
get. 


Bi  = 


_ dx _ 


.  . .  (II- 8) 


and  finally, 


K  fluid 
K  solid 


s 


^  _  K  solid  6  Re’^Pr 
^  ^  K  fluid  r  3":3r(4>i 


.  .  .  (U-9) 

.  .  .(II- 10) 


r  ^ _ djc _ 


...(II-ll) 


The  use  of  eq.  (2.  14)  instead  of  (2.  12)  in  the  development  above  is  unfortunately 

not  possible  as  it  leads  to  an  equation  which  is  not  easily  soluble.  Solution  of 
(*) 

(II- 1 1 )  therefore  proceeds  best  numerically.  This  was  carried  out  on  the  Stanford 
B5000  digital  computer  .  The  program  makes  use  of  the  Fourth  Order  Adams - 

^  ^  Using  (2.  12)  to  evaluate  the  integral. 


1 
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r  1  (*) 

Bashforth  Predictor  -  Corrector  Kif-thod  for  the  numerical  integration  [3 3j  ,  , 

and  results  are  given  below  and  in  figure  (19),  and  table  (I). 

The  equation  to  be  solved  is, 

T"  -  +1=0  ...  (II- 12) 

X 


with, 
T  = 


e  = 


^fluid 

solid 


solid 

^fluid 


© 

M 


...(II- 13) 


...(11-14) 


and  boundary  conditions, 
T'(0)  =  T'(l)  =  0 


...(11-15) 


As  explained  above,  the  importance  of  the  second  derivative  term  in  the  equation 
(II- 12)  will  be  felt  in  narrow  "boundary  layers"  near  x  =  0  and  x  =  1.  Near  these 
boundary  points  all  terms  in  the  equation  should  become  of  comparable  order  of 
magnitude.  This  may  be  emphasized  by  introducing  streched  cobrdinates  as 
follows: 


A 

T 


=  £*^T: 


=  6  x; 


...(11-16) 


The  boundary  x  =  0  coincides  with  x  =  0  for  any  non-zero  £  .  V/ e  thus  see 

i/s(  +  ) 

immediately  that  near  that  boundary  T«>c  6  .  Table  II  below  gives  some 


numerical  data: 


TABLE  II 


numerical  integrtn. 

{o.zo/e  y  ^ 


] 


0.  20 
1. 00000 
1. 00000 


0.  15 
1. 03737 
1. 0391 


0.  10 

I.  09867 

J .  1488 


0.  05 
1. 21640 
1.  3195 


0.001 
1.  9525^ 
1.  810 


/  i^\ 

'  The  author  is  indebted  to  Mr.  Ury  Passy  for  help  with  the  numerical  analysis. 
This  was  pointed  out  to  the  author  by  t)r.  A.  Acrivos. 
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PROGRAM  FOR  NUMERICAL  INTEGRATION 


SIANFURD  dSOCO  ALGOL  -*  VERSION  216/6a 

’lEGlN  CUHMERI  ROrtR  2EEV  COMIuCTlQN  E^FLCI  CALClLATIQNS  FOR  THE  CONSTANT 

Flu*  CASL  FUK  various  HIUT  number  values  bY  AUAmS  (jASRFQHTH  integhaticn; 
REAL  2LrA,EPSllINIEGER  K»lJHE4L  ARRAY  APR . *F L I  1 1 2  1  *  I M T 1 AL 1 » F 1 N AL2 
lOtlOOJ/LABEL  ENUltEN02>UUAGAiN; 

►  QRHAT  FHI 1(///»"IHE  VALUE  OF  EPSlLCN  IS  NUn " » X 2 » f 5 . 3 )  I 
PROCEDURE  FUNCT(NU.XPH,F);v4LuE  NUFREAL  NUiHEAL  ARRAY  XPR[*J»FC*j; 

BEGIN  FtlJ«.XPR[2JFF[2U(XPRCl]/NU*(l/i)-n/EPSi;END  FlNCn 
REAL  array  KFURAUAMSIO  I  3»  1 1  30  J  »  Y  INCF'jR  A04MSC 1 1  30  1 J  CC“PENT  GOES 
BEFURE  AUAHS3 

procedure  AUAPSCSI2E»  w2ERO»  INITIAL.  FINAL.  PRIM.  RELH»  ABSd. 

YIMTIAL,  YFINAL*  FUNCT)  i 

VALUE  SUE.  MZEHC.  INITIAL#  FINAL.  PRIM.  RELri.  APSB  > 

INTEGER  SUE  I  REAL  H2ERC.  INITIAL.  FINAL*  PRIM.  RELB.  ABS0  ) 
real  array  YISITIAL.  YFInAHU  > 

PROCEDURE  FUNCT  I 

begin  COHHENT  ADAMS  VERSION  OF  APRIL  1.  1V6A  I 


RETUMM  END  ADAMS  ) 

►  UR  £PSI»5R-2  STEP  5P-2  UNTIL  0.3  00  BEGIN  I  «-0 1 nR  I T E ( FM  T 1 .  t P S I  > I 
xPR[ll«'5*i‘-llXFRI2l*CI 

IN|IIALUIl*XPfi(l]IK*Ol2ErA*0.llFOR  Ul  STEP  I  UMIL  lOO  00  BEGIN 
UOAGAJNI  ADAMS(2.0.01.P-l0.l.2ETA.?-?.P-2.XPR.XFL.FUNCI)l 
IF  ABS(xFL[2J)>  I'-6  THEN  UEGIn  IF  I»1  THEN  BEGIN  F  I NAL2 1  i  •  1  I  y-xF  L  t  2  I ) 
XPRl  1  INITIALK  I  JMMTIALtC  I-n*?-/j;GO  TO  ENOllENO  ELSE  BEGIN 

►  lNAL2t  l-n*XFL[2JIXFR(  1  Jy-INITIAL  in  IMNII  lALlC  I-l  ]-(  INITIAL!  t  I-l  1 
-:NIIlALUI-2J)/<FINAL2(I-lJTlNAL2tI-2  1)>«FINAL2CI-niGG  10  ENDIHNOI 
ENOMF  K*0  IHEN  BEGIN  K«.  I  i  ZET  ay-0 . 05 1  GO  TO  OOAGAIN 'ENOMF  KAO  IPEN  GO  TO 
EN02IEND1 lEN0ICN02tEN0IEN0. 


GRAPH  FIGURE 
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Near  the  other  boundary  the  investigation  of  the  properties  of  the  so¬ 
lution  proceeds  as  follows;  the  complementary  solution  to  (II- 12)  is, 


T  = 


t  6  a5/6,  ^  „ 


.  6  -*5/6.1 

3/5'-r’‘ 


.(II-H) 


with'X  =  X.  £  -  [the  value  of  the  particular  solution  of  the  full  equation  is 

approximately  unity  at  x  =  1  (i.e.  x-*kx)  ).  Now,  as  the  function  diverges 

as  the  argument  increases  without  limit  we  must  have  c  =  0;  also,  an  asymp- 

(*)  ^ 

totic  expansion  may  be  used  for  the  function  K_,_.  The  result  is  then  easily 

3/5 

derived  that 


1  -  Tor.e^/^.  exp(-  1  iji) 


. .  .(n-i8) 


This  agrees  qualitatively  with  the  results  of  the  numerical  integration. 

1  / 

Standard  singular  perturbation  analysis  would,  of  course,  yield  ATatE 
using  the  standard  stretching  technique.  From  the  numerical  integration  it 
appears  however  that  for  identical  values  of  8  ratios,  the  ratio  AT,  /  A 
depends  on  the  absolute  value  of  also;  of  this  the  standard  analysis  gives 
no  indication. 


TABLE  (III).  Numerical  Solution  Near  x=l 


€ 

T 

1  -  T 

0.  05 

0. 911074 

0.  088926 

1. 000  -j 

0.  10 

0. 869471 

0.  130529 

,  h  1.468 

1 . 468  J- 

0.  15 

0. 841093 

0.  158907 

i.  788  >■  1.  374 

0.  20 

0. 82065 

0. 17935 

2. 016  J 

(«) 


See  c ,  f.  [34j ,  p.  1 50. 
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It  is  also  possible  to  develop  a  simple  integral  approach  for  estimating 
the  role  of  circumferential  conduction. 

Introduce  the  assumption  of  constant  mean  ^(x)  into  equation  (2.  11);  then, 


_  ,  .  ,,  r  Pr'^Re-t.  ,  ..^3 


whence 


■19s(x.) _ !_ 


(3y3y.  rci) 


.  (U-  19  ) 


.  .  (II-2Q) 


The  circumferential  flux  per  unit  axial  length  of  cylinder  is  then, 


-K  solid 


£  d9g(x  )  _  solid  Pr  Ret. 


K  fluid  cs^j^rd)' " 


'^3 


...(n-21) 


The  radial  flux  per  unit  axial  length  of  cylinder,  for  a  strip  of  width  (x  -x*!')  is. 


q"  r  (xj_-x’!') 


.  .  .(11-22) 


Thus,  a  measure  of  the  relative  magnitude  of  circumferential  to  axial  fluxes 
is  obtained  from  (11-14)  and  (11-15); 


1  _ 1 

m  (Sgyi 


C  ^ 

K  solid  J>_  Pr  Re  _ 

K  fluid  r  (x^-x*)  *43 


.  .  .(11-23) 


i.e.  the  reciprocal  of  the  Biot  number  is  a  direct  estimate  of  the  importance 


of  circumferei-tial  flux  compared  to  radial  flux. 


Finally,  equation  (II- 11)  gives  the  possibility  of  estimating  the  departure  of 
an  almost  constant-flux  surface  from  the  ideal  case  when  the  parameter  £  anc. 
the  measured  values  of  temperature  are  known. 

Rewriting  eq.  (II- 11)  thus. 


er-  >  =  0 


.  .  .(11-24) 


E«l 


leads  to. 


^  i  {:3c)  +E  j  (y:)  (x)  +  ...(11-25) 

where  £iJC)  is  obtained  from  the  measured  variation  of  HP.  and  is  then 

obtained  through  numerical  differentiation  of  these  values.  Provided  a  sufficient 
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number  of  data  points  are  taken,  a  first  order  correction  can  therefore  be 
obtained.  This  is,  however,  liable  to  be  a  rather  inaccurate  procedure.  Now, 
it  is  easily  shown  that, 


Nu  = 


.  .  .(11-26) 


therefore 

Nu 


r  I  ciillxj 
Ucj:)  ^  J-Cx) 


.  .  .(11-27) 
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Appendix  III 


Phygical  constants  of  electrolytic  solution 

There  are  wide  discrepancies  behA/een  values  of  the  physical  prop¬ 
erties  of  our  electrolytic  test  solution,  as  found  b/  various  investigators. 
We  were  therefore  prompted  to  compare  these  values  in  figures  (21)  through 
(24). 

It  appears  from  an  examination  of  figure  (24)  that  both  the  composition 
and  concentration  of  supporting  electrolyte  seem  to  have  a  notable  influence 
on  D,  as  well  as  the  manner  of  testing. 

The  one  systematice  trend  discernible  seems  to  be  the  decrease  in  D 
as  supporting  electrolyte  is  added.  An  unexpected  feature  is  the  fact  that 
'^(K  Fe(oN)  )/0{K  Fe(CrN)  )  can  be  either  largeror  smaller  than  unity, 
aepending  on  the  temperature  and  on  the  concentration  of  supporting  electro¬ 
lyte.  It  also  appears  from  ^35j  that  diffusion  coefficients  found  using  KCl 
as  a  supporting  electrolyte  are  somewhat  higher  than  those  measured 
with  NaOH  as  supporting  electrolyte. 


Appendix  IV 

Data  reduction  of  heat  transfer  experiments 

For  data  reduction  the  Algol  p»‘Ogram  of  ref.  |^36]was  used:  this  program 
includes  a  correction  taking  into  av  i  ourit  the  dependence  of  fluid  viscosity  c  i 
temperature  based  on  [bt].  The  sole  modification  introduced  into  the  program 
was  a  correction  for  the  dependence  of  the  electrical  resistivity  of  the  copper 
heading  atrip  on  temperature.  A  printout  of  the  program  is  given  on  page  60. 


i  (0 


i\o'J  Xbhi'ifAtJ’t 


DIFFU 
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»IUlt03*  1, rsoo-c 9, <75796 1-<  6, 0d50639*C  5. 229928 1-{  5665. 288-(  2*5, 5829ft 2 
•3,79ft65«9>«THt01  )KTHr01)*THC3lJCl,08-ft;««THtOJ  )KTHt 03 *1 .0#-2 ) 

>(THC03Kl.08*2i 

f«St03*2ftl.90«8KC7t33*MUSt03/KTHC33J 

PRS3t0]*7A5t03*(l/3}} 

fACT0lll03«‘NUt03/>»RS3t03i 

IfACTOI^t  03*(  1-(1 -Must  03/MJtj3  >'<0.095576  )xr*CT09C  03) 

CNCi 


RRITt(rMTl2frOR 

0«-l 

STEP 

1 

UNTIL 

P 

00 

TSt03)l 

llRITE(rMTl3#rOR 

3<-l 

STEP 

1 

UNTIL 

p 

DO 

CPt03)l 

NRlTC(rMTlft*rOR 

0«-l 

STEP 

1 

UNTIL 

P 

00 

KTHt03 )) 

ilRlTC(rMTl5«rOR 

3*>1 

STEP 

1 

UNTIL 

P 

00 

NUt03 )) 

•IRIT£(FMTl6,rOR 

0*1 

STEP 

1 

UNTIL 

P 

00 

PRSt03>) 

i(RITE(FMTl7»rOR 

Q*l 

STEP 

1 

UNTIL 

P 

DO 

PRS3t03>> 

»<RITE(FMT18»rOR 

0*l 

step 

1 

until 

» 

00 

rACT0Rt03>> 

i<mTC(fMTl9#rOR 

0*-l 

STEP 

1 

UNTIL 

P 

DO 

GFACTORtOI)! 

HRIT£(rMT24)j 

e«iOJ 

HRlT|(rMT23»K)l 
SO  TO  LEOril 

IHRlTECrMT9)  j  GO  TC  LE0r;MlKAN2l>lRITECrMTl0>>G0  TO 
^XKAM3lWRITCCrMT20))  GO  TO  |.EOri  MlKANft l  i»<RlTE{rMT22) I  SO  TO  LECF3 
UCOrilENOiENOi 
ucoriENO. 
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»J38  T221  ALSOl  2  1200  '<0TE«»09  ZEEV  293  Nj  ST«1P  CYLIN0E9 

begin  comment  ROTEH  ZEEV  calculation  op  NU  pop  strip  CTlINOERJ 

INTEGER  1»  J»‘<»Ll»L2»l.3#LA»‘'»N,0,P,a) 

INTEGER  ARRAY  DMPfiA 1 1 1 50 1 ; 

real  array  v»  n  J»  TR  JLKa  JONE  »  R0»  RPYN»  RE-^  alP  »  VOlT  age#  ThpT  a  1 1  1 50 1 » 
TS»KTH»CP’»NU»PRS.RRS3»’^US.PACT0R»GPACT0R»Ti^,RESC0Rtll200]) 
label  STARTaNMANlaMlKANZ.wnANOaMT^  ANA>LE^''^LE3P1  J 

PORHAT  rNTH/Z/aKlSa-RESULTS  QP  RUN" »  X?  ,  I  A  .  X  3,  -  ARE  AS  FOLlOwS"# // ) ) 
format  rNTZl^O^EjA*"*  XA» IPIO, 3)»PmT3("T3UL<*"» YA» iFlO, 3)» 
rMTa{"Ul*''»x7»iri0,3)»PMT5C’R0«*’»X9»lF9.A)» 
rMT6{  POISES)  ■"»  X3»  ire.  5  )»PMT7("RE«"»X7»ri0.3)» 
rNT8(-REMALP«"»X3,iri0.3)»  ^ 

PMT9{  "NUMBER  OP  3aTA  POINTS  ilP  TBUl^  Y  NUmRER  OP  DATA  POINTS  QP  3MEGA  )» 
PmT10(**NJMBER  data  points  volts  *  THAT  POR  OMEGA"), 
PMTU(*’V0LTAGE«*’»X6,ePfl.5), 

rMTl2("TS«**,X8,8P9.3,X2),PMT13("CP»’’»xe,8P9.3,x2), 

rMTH("<«*’»Xl0,8PiJ.A,Xl),'*^Tl5("NU»"»Y8»8P9.3»X2), 

rM?l6("PRS«",X7,9r9.3»X2)»PMTl7("PRl/3»",X5»SP9,3,X2), 

PMTlB("NU/(PRl/3)«",8P9,3,X2),PMTl9("3NJ/PRl/3«",Xl,8P9,3,X2,//), 

rMT20("N3  CARO  POR  m  OR  Ll>20  OR  l2>20  OR  NO  NJMB  SmITCM  »0S"), 

PMT2H  "Ll*"»  IR )» ^‘^^221  "O^EOA  OR  TRJLM  0AT'J'<  MAS  VALUE  ZERO")* 

PMT23( /,X15»  "END  OP  R UN**,  XI ,  I A  ) ,  PmT2R  ( // )  1 
STARTIBEGIN  read  (M); 

IP  M.o  Then  go  to  MIXAN3;  ip  m  >12  then  go  to  mi<an3; 

POR  N*-l  STEP  I  JNTI^  DO  BEGIN 
REAO(<»Ll»POR  1“!  STEP  I  JnTIl  LI  00  O'^EGaCID) 

RE*0(L2»P0R  I*l  STEP  I  until  L2  00  TBJLMtlDl 

RE*D{l3»P0R  I*-l  STEP  I  until  l3  00  VOlTAGEIIDT 

IP  Ll  *  1.2  Then  go  to  mIkanUTp  lT>23  ^^EN  GO  TO  mi<an3; 

POR  I*-l  STEP  1  until  Ll  00  BEGIN 


IP  L3  *  Ll  then  go  to  mikanz; 

IP  0MEGA[I3«0  Then  GO  TO  m!<».nA;IP  TgjLMm»0  then  go 
HRITEIP^TI.O;  HRITE(PMT21,LI  )1 


TO  mIKANAIEND; 


POR  J*l  STPP  I  until  Ll  '^0  BEGIN 

YtJ3*-(3MEGA[J]*3*P-6)/25.T73I 


THETAI  J3*>TBJLXt  J3-2O.OO; 

MUt J3*l.Z800-(9,A75Z96l-{6,0650639*(5.229929l-{5665.288-(2A5.5929A2 
•3,79A65A9xTHETA[ J3 )»THETAC J3 IxThETAC J3*l ,OP-C)*THETAf  J3 )* 
THETAC  J3 -‘1.08-2  I-TMET  A  [J3  *1,08-21 
UONEt  J3*-12.25>‘S3RT(v[  j3/MJt  J3  )J 

ROt  J3*'1.0R8330-{I0.030225*T3UL'<C  J3*A.  306)>‘l.0?-fi)>‘T9UL'<t  JIT 

REYNI  J3*-l  .ROSwJONEI  JT^ROI  JT7MUC  J3 1 


REHALPI J3*SQRT{REYN[ J3 )T  ^ 

HRiTe(PMT2»0MEGA[  J3  )JHRITE(PMT3,T8JL<C  J3)T  MRITE(PMTA,uONEt  JDT 

HRlTE(PMT5,R3[J3);HRITE(P“T8,MU[J3)TMRITE(PMT7,REYNtj))) 

HR 1TE(PMT8»REHALP[J3)1HRTTE(PMT11, VOLTA GEtJ))T 

READ(P»P0R  a*-l  STEP  1  UNTIL  P  00  TSCQ3)I 
IP  p>20  Then  go  to  mixans; 

POR  0*-l  STEP  I  UNTIL  P  DO  BEGIN 

CPt03*0,A30*S.D098-A)‘CT5[03-l0)I 

XTH[03*(9,262-68-Ak{TSC03-10))*1.0*-21 


RE5C0Rt03  *'l*{3.9  3?-3)*{TSt03-23.1  )I 

NUt01*-6,69l63572»AwV0LTASE[J3*2/<(TS[03-T3UL'<tJl)“RESC0Rt93)J 


THt03*-TS[03-20.001 


